Deuteron Matrix Elements in Chiral Effective Theory at Leading Order 
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We consider matrix elements of two-nucleon operators that arise in chiral effective theories of the 
two-nucleon system. Generically, the short-distance piece of these operators scales as 1 /r™ , with r 
the relative separation of the two nucleons. We show that, when evaluated between the leading- 
order wave functions obtained in this effective theory, these two-nucleon operators are independent 
of the cutoff used to renormalize the two-body problem for n — 1 and 2. However, for n > 3 general 
arguments about the short-distance behavior of the leading-order deuteron wave function show that 
the matrix element will diverge. 



Introduction: In recent years chiral effective theory 
(%ET) has achieved significant prominence as a technique 
via which model-independent results can be obtained 
for few-nucleon systems (see Refs. 0, for reviews). 
XET is based on the realization that the interaction be- 
tween pions and nucleons is governed by the approximate 
chiral symmetry of QCD. Therefore the use of heavy- 
baryon chiral perturbation theory, together with Wein- 
berg's proposal that in few-nucleon systems the quantity 
that has a well-behaved chiral expansion is the nucleon- 
nucleon potential V 0-13. seems to facilitate a systematic 
calculation with which nuclcon-nuclcon scattering can 
be well described — as long as the collision energies are 
significantly below the chiral-symmetry breaking scale 
A x [sj. liJ. Lj] In this approach V — and also by extension 
two-nucleon-irreducible operators for other processes — is 
expanded as a chiral scries in the usual chiral perturba- 
tion theory expansion parameter 



P : 



p, m-x 
A X ,M' 



(1) 



where is the pion mass, p the momenta of the nucle- 
ons involved in the scattering, and M the nucleon mass. 
Although questions have been raised about the consis- 
tency of such a power counting [I, El, El El this approach 
has had considerable success in describing the scattering 
data in the two- and three-nucleon sector. 

The theory has also been shown to be consistent — 
in the renormalization sense — in the 3 Si- 3 Di channel, 
the channel where the two-nucleon bound state deu- 
terium occurs 0, Ih3 ] . The leading-order wave function 
of deuterium, \ipho) can therefore be obtained by solving 
the Schrodinger equation — in either momentum-space or 
configuration-space — for two nucleons interacting via the 
piece of the NN potential which is of chiral order zero. 
As first realized by Weinberg 0, that piece consists of 
the venerable one-pion-exchange potential, together with 
a (momentum-independent) four-nucleon contact inter- 
action. This potential is singular: the Hamiltonian it 
generates is unbounded from below. It therefore requires 



regularization and renormalization. In practice, the po- 
tential is regulated at some momentum scale A (or in co- 
ordinate space at a distance 1/A). The strength of the 
contact interaction is then adjusted to reproduce some 
observable, usually the deuteron binding energy. If this 
can be done over a wide range of A, and if other NN scat- 
tering observables are independent of A up to corrections 
of higher order in the x^T, then we say that the poten- 
tial has been renormalized. In Ref. Beane and collab- 
orators showed that this could be done for the 3 Si- 3 Di 
channel — at least if one was only concerned about the 
results in the chiral limit m n — 0. This conclusion has 
since also been reached using a different regulator and in 
a regulator- independent fashion in Refs. [llj, liJ] • 

With these wave functions in hand it is natural to re- 
examine the many successes that phenomenological po- 
tential models have had in describing the deuteron's in- 
teraction with external probes, such as electrons, pions, 
or photons. In such a calculation we use an operator O 
that is appropriate to the particular external probe un- 
der consideration, and has been derived using a chiral 
expansion in powers of P. That operator is then sand- 
wiched between a wave function obtained from a chiral 
potential (which presumably should be computed to the 
same relative order), yielding matrix elements: 



M 



101 



(2) 



Many such calculations have been performed (see 
Ref. [lij ] for a recent, partial summary), although, with 
one exception which we will discuss below, none of them 
have employed the wave functions of Refs. 0,0,0]- In 
these calculations the bound-state wave function is found 
by employing some regularization prescription involving 
a cutoff A. Here we will discuss the conditions under 
which the use of such wave functions in the evaluation 
of Eq. yields a matrix element that is independent of 
A. The extent of the A dependence in the result for M. 
tells us the degree to which the %ET prediction for the 
matrix element is model independent. 

Recently, Meifiner et al. [l£j and Nogga and Hanhart 



[lif considered this question for pion-deuteron scatter- 
ing, a reaction originally discussed within this frame- 
work by Weinberg in 1992 17]. There a two-nucleon 
operator representing the process depicted in Fig. ^a) 
gives a large contribution to a n d- The presence of this 
large double-scattering term can obscure attempts to ex- 
tract the isoscalar pion-nucleon scattering length from 
a^d- Pion scattering on the individual nucleons in the 
deuterium nucleus (see Fig. ^b)) is 0(P 2 ) in the yET, 
and yields a piece of the 7rd scattering length 0> uM '■ 



i (b} 



(i + aO 

(l + M/2) 



(3) 



with /i = m v /M. Meanwhile the double-scattering dia- 
gram of Fig. [Ha) gives [H El El : 



,(«) 

l Ttd 



4^(1 + M /2) V2/2 




where here, and throughout what follows, 



</(r)> 



dr/(r)(u 2 (r) + w 2 (r)), 



(4) 



(5) 



with m and it; the 3 Si and 3 Di radial deuteron wave func- 
tions. 

The contribution of Eq. Q to the pion-deuteron scat- 
tering length is suppressed by one power of P relative to 
the nominally- leading contribution @. However, since 
corrections to the leading-order NN potential are sup- 
pressed by two powers of P 0, a next-to- leading order 
computation of a^d involves the evaluation of the matrix 
element: 



(6) 



LO,A 



where the subscript LO indicates that the matrix ele- 
ment should be taken with deuteron wave functions ob- 
tained with the leading-order xPT NN potential, and 
the subscript A is included in Eq. JBJ) because a cutoff 
must be applied to that potential before it is renormal- 
ized to give the correct deuteron binding energy. Since 
there are no short-distance pieces of the 7r-d operator O 
at next-to-leading order the matrix element © must be 
cutoff independent over a significant A range if Wein- 
berg's approach is to be a consistent way to calculate 
pion-deuteron scattering. 

Nogga and Hanhart gave numerical evidence that the 
matrix element (0 has a A — + oo limit that exists and is 
finite. (This is in accord with a similar conclusion previ- 
ously obtained in Ref. where (-) was evaluated us- 
ing wave functions derived solely from one-pion exchange. 
The details of the evaluation of Ref. E3 will be explained 
further below.) This shows, Nogga and Hanhart claim, 
that expanding O in powers of P is a consistent way 



(a) 



(6) 



(c) 



FIG. 1: Three Feynman diagrams for different contributions 
to the 7rd scattering length in chiral effective theory. The dots 



are vertices from C^, while the sliced blob is a vertex from 



,(2) 



to calculate the pion-deuteron scattering length. Con- 
sequently, at least in principle, it should be possible to 
perform a high-accuracy, model-independent extraction 
of the pion-neutron scattering length from pion-deuteron 
scattering data. 

But other two-nucleon operators that occur in 7rd scat- 
tering computations have higher powers of r in the de- 
nominator in their co-ordinate space form. For instance, 
the next term in the 7rN multiple-scattering series, de- 
picted in Fig. Ofc) , yields a contribution to a^d [H3 : 



Ac) 



1 



16tt3(1 + M /2) V2/2 



l 



(7) 



This contribution is of order P 5 , and so is N 3 LO in the 
chiral expansion for a^d, but numerically it is the next 
two-nucleon effect that must be considered after the in- 
clusion of the matrix element iJBJ 20] . The failure of the 
standard chiral expansion to account for the large size of 
such effects is related to the absence of any suppression 
of the operator JJJ at long distances. It led the authors 
of Ref. 20] to propose a different power counting for the 
7rd scattering length — one that more accurately captures 
the relative hierarchy of mechanisms contributing to a I( j. 
Regardless of what the correct counting for the operator 
l [7| l is, if the impact of the triple-scattering mechanism 
of Fig. njc) on extractions of the isoscalar pion-nucleon 
scattering length from a^d is to be assessed we must not 
only consider the matrix element JSJ, but also: 



l 



(8) 



LO,A 



More generally, all two-nucleon operators in the chiral 
effective theory will scale as 1/r" at distances <C l/m^ 
and C 1/? (with q the magnitude of any momentum or 
energy transferred to the nucleus by the probe). The 
value of n is determined by the process under consid- 
eration and the order to which the two-body operator 
is calculated, with higher values of n being reached as 
higher orders are computed in the chiral expansion for 
O. In a perturbative calculation of deuteron matrix ele- 
ments within xET the highest-order (and most singular) 
pieces of the operator will contribute to the overall result 



3 



for M. through their expectation value taken with the 
leading-order wave function. Thus, we are led to a very 
general question: If one considers the expression 



(9) 



LO,A 



then for which values of n does the A — ► oo limit exist, 
and for which is it finite? It is this question, together 
with associated ones involving operators that connect 
the S- and D-wave components of the deuteron wave 
function, that we will answer in this paper. 

Theory: Nogga and Hanhart solved the Schrodinger 
equation in its momentum-space form, i.e. the homoge- 
neous Lippmann-Schwinger equation: 

(p|^lo)a = G (p) / A 7 ^3^ (0) (p,p')<p'IV'lo)a , 
Jo ( Z7r ) 

(10) 

where A is the scale at which the potential V is regulated, 
and Gq(p) — (—Bd—p 2 /M)~ 1 is the (free, center-of-mass 
frame) two-nucleon propagator, with B d and M denoting 
the deuteron binding energy and nucleon mass, respec- 
tively. The leading-order potential is given by a one-pion 
exchange (OPE) contribution and a short-distance piece: 



V(°'(q) 



9A_ 

2f« 



n • r 2 - 



(ax ■ q)(cr 2 • q) 



-in 



C t P t , 
(11) 



with q = p' — p the three-momentum of the exchanged 
pion. In Eq. P t is a projection operator that projects 
onto the 3 Si channel, and Ct the strength of the short- 
distance potential in that channel. While the OPE con- 
tribution is totally determined at leading order through 
the pion mass m n , the axial coupling constant g A = 1.26, 
and the pion-decay constant f„ — 92.4 MeV, the contact 
interaction parameter Ct has to be determined from NN 
data and will be a function of the cutoff A. 

From Eqs. I|1U|) and l|ll|) it is straightforward to obtain 
the coupled one-dimensional differential equations which 
describe the deuteron wave function: 

-u"(r) + U s (r)u(r) + U sd (r)w(r) = -~/ 2 u(r), 



-w"(r) + U sd (r)u(r) + 



w(r) — — 7 2 w(r), 



where u and w are, as defined above, the deuteron radial 
wave functions. The coupled-channel potential is given 
by 

U S = U C , U sd = 2V2U T , U d = U c ~2U T , (13) 



with 

U c = 
U T = 



m\Mg\ e~ m - 
16tt/2 r 

m\Mg\ e'™-- 
16tt/2 r 



m^r (m T r)' 



(14) 



Equations fT^ - lfTl^ will be equivalent to Eqs. ((THf and 
(tTT|l provided that r > 1/A. The equations (fT2*)l are 
solved subject to the following boundary conditions as 
r — > oo: 



A 



-77" 



u(r) 

w(r) -> riA s e~ ir ( 1 + — + 



7r ("fry 



(15) 



with 7 = \JMB d the deuteron wave number, As the 
normalization constant which guarantees that 



dr(u 2 (r) + w 2 (r)) = 1 



(16) 



and rj the asymptotic D/S ratio. 

By employing Eqs. (|12fl 114(1 . and solving them for ar- 
bitrarily short distances r, Pavon Valderrama and Ruiz 
Arriola have calculated u and w using boundary condi- 
tions at inter-nucleon distances of order 0.1-0.2 fm [l2j |. 
This echoes the much earlier work of Sprung and collab- 
orators [U, as well as the solution of the leading-order 
potential with a square- well regulating the short-distance 
behavior [£j . Within such an approach, an analysis of the 
asymptotic short-distance behavior of the components u 
and w gives the following result |J 0, : 



u sd (r) = A, 



1 (r 



3/4 



w sd (r) 



A, 



1 / r 

V5\R 



+2 3 / 2 |C 2A |cos(4^/- + ( / ) 

3/4 



+2\C 2A \cos 4W- 



(17) 



C 2 a and C2R. are normalization constants which have 
been determined in 01 j is a new length scale that 
enters the non-perturbative problem. It is defined by 



(12) R 



3g 2 A M 



When Eqs. JEJ-djl are solved in Ref. [13] 
> is determined by the boundary condition 
and so (j) 1S regulator independent, and is 

7, and R. After 



the phase 
at r = 

a function only of the scales m^, 



computing the numerical solution of Eqs. I(12|l to a 
sufficiently small radius r, it can be matched to the 
r — > form of the deuteron wave function l|17|) and 



4 



an — in principle — regulator-independent wave function 
can be obtained 1 . 

Results: Let us now return to the problem of com- 
puting 



(18) 




1 



First, we split the integral up as: 

* u 2 ir) + w 2 (r) 

dr 

o r n 

R * ^ u 2 (r)+w 2 {r) r°° u 2 {r)+w 2 (r) 

T r 

(19) 

where R* is sufficiently small that the asymptotic forms 
(|17|) apply, but is still ^> 1/A. The piece of the integral 
from i?* to oo can be calculated within xET and will 
depend only on low-energy scales such as 7 and m w , and, 
of course, on R* itself. Meanwhile, as A — > 00 the first 
integral goes to zero, as long as the integrand is inte- 
grable. Note that any regulator dependence of the wave 
functions should be contained in that first piece of our 
integral, since we have already explained that the wave 
functions for r > 1/A are regulator- independent solu- 
tions of Eqs. <|12jl - 1114 (1 . Therefore, in order to establish 
whether or not (l/r n )\ is regulator-independent the per- 
tinent piece of Eq. 119J1 is 



1 



sd 



dr 



u 2 (r) + w 2 (r) 



(20) 



If this has a A — > 00 limit that is finite then the entire 
matrix element (l/r")^ will also be well behaved in that 
limit. Inside the integral Eq. ljUJi we may substitute the 
expressions Eq. l|17l) for u and w. We then see that short- 
distance integrals involving the exponential piece of the 
wave function will always converge for any n, as the ex- 
ponential itself regularizes the result. On the other hand, 
the term including the cosine function does not vanish at 
the origin and therefore we have to determine the values 
of n for which the integral 



1! 



drr 3/2 ~ n < 



1/A 




(21) 



gives a finite result as A — > 00. 

This question is easily answered by simple dimensional 
analysis. For n = 3 or higher the integral in Eq. 121|) 
will diverge. However, for n = 1, 2 if we ignore the ex- 
ponential pieces of u and w we can solve the integrals 
analytically and obtain 




where: 



f n (x) = / dtt 2n ~ b cos(t + 2</>) 



(23) 



The function f n can be written as a linear combination of 
incomplete Gamma functions. Its asymptotic expansion 
is: 

f n {x) - a n + h n (a)+x 2n - 6 sm(x+2<f>) + 0{x 2n - 7 ), (24) 
cis x - — > oo (n = 1, 2), with: 

7T 7T 

= — cos(20); a 2 = -- cos(20). (25) 

The function h n (a) can easily be evaluated, but its value 
is not relevant for our purposes here, since it cancels in 
Eq. I|22|) . Therefore for n — 1 and 2 the matrix element 
has a finite A — > 00 limit. However, that limit is ap- 
proached in the presence of cutoff-dependent oscillations 
whose amplitude scales as (i?A)™~ 3 . Another cycle in this 
oscillation occurs whenever 1/A becomes small enough 
that we integrate over another node in the deuteron wave 
function. The presence of oscillations is thus associated 
with the appearance of spurious bound states in the ef- 
fective theory — bound states whose binding energies are 
greater than the theory's breakdown scale. However, for 
the n = 1 and n — 2 matrix elements, the oscillations 
vanish as A — > 00, so even though such bound states are 
present, the procedure of evaluating the matrix element 
with a wave function at a given A should yield the correct 
(i.e. A — > 00) answer, as long as a sufficiently high cutoff 
is chosen. 

Now adding the A — > 00 piece of these matrix elements 
to the piece that comes from integration between R* and 
infinity, we find: 



1 Although the final wave function is in principle regulator in- 
dependent, the boundary conditions employed in the numerical 
computation of the long-range components can lead to cutoff ef- 
fects at short distances and this leads to a small numerical uncer- 
tainty of the constants C2A, C2R quoted in |l2j. However, with 
sufficient care this uncertainty can be made arbitrarily small. 




= 0.478 fm" 1 , 



0.425 fm -2 , 



(26) 
(27) 
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FIG. 2: Cutoff-dependence of the matrix element (Ma, using 
the radial deuteron wave functions employed in The 
cross is the result found using the wave function of Ref. [l2^ . 
as given in Eq. (1261 . 
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FIG. 3: Cutoff-dependence of ( jr)A, using the radial deuteron 
wave functions employed in [lrl | . The cross is the result found 
using the wave function of Ref. as given in Eq. 1271 . 



where the first result is in agreement with the number 
given in Ref. [1^, As pointed out there, it is essen- 
tially consistent with the range 0.450-0.465 obtained for 
(I/r) us ing a variety of potential-model wave functions 
in Ref. |20j. The result l|27|l is new, and does not fall 
within the range 0.286-0.345 fm 2 quoted in Ref. [Ij. It 
takes longer to reach asymptopia for this matrix element. 

In order to check the results of Eqs. (|26|l and l|27|) we 



have computed the expectation value (V'l 



-IV'lo) for 

n = 1 and 2. We used the same wave functions as were 
used in 0] with cutoff values between 2 and 20 fm -1 . 
As Fig.[2]shows, we find similar behavior for (1/r) as was 
displayed in that paper, and a limiting value that is con- 
sistent with that found by Nogga and Hanhart [22]. But, 
we can now interpret the A-dependent oscillations seen in 
Fig. [21 as exactly the ones predicted by Eq. I|22|l . Mean- 
while, Fig. |21 shows that, as expected based on our anal- 
ysis of Eq. (|22|l . these oscillations are more pronounced 
for (1/r 2 ), and in consequence a much higher cutoff is 
needed to achieve a converged result. However, also for 
this matrix element the numerical results indicate that a 
finite limiting value exists, which is in accord with our 
previously presented analytic arguments. 

Since Eq. (|17|) identifies the relevant behavior of the 
r — * piece of the deuteron wave function, we are able 
to predict that all powers n > 3 will lead to divergent re- 
sults. This result is supported by numerical calculations 
of the corresponding matrix elements with the wave func- 
tions used in [l6| . 

We also find that, to a very good approximation, the 
result for 




(28) 



LO.A 



agrees (for n = 1 and 2) with the numerical results 
obtained with the co-ordinate space wave function of 



Ref. [13 that we presented in Eqs. (2J) and (The 
slight disagreement exhibited in Figs. El and 01 between 
the trend of the dots and the cross that represents the 
result of Eqs. (|26|l and l|27|) can presumably be traced to 
the 3% difference in the 7rNN coupling constants used in 
Refs. [12 and 0.) This is an important result, because, 
together with the even better agreement for matrix ele- 
ments of positive powers of r (see Table I of Refs. 0,0) 
it suggests that the procedure of constructing a limiting 
sequence of cutoff-dependent wave functions and evaluat- 
ing the matrix element as a function of A is not actually 
necessary. Instead such matrix elements can be com- 
puted using the co-ordinate space results for the deuteron 
wave function given in Refs. dEHHl. 
Next we consider integrals of the form: 



dr 



u(r)w(r) 



(29) 



For n — 1 this integral occurs in other, numerically less 
important, corrections to the 7rd scattering length. For 
n = 2 contributions of the form (|29|l occur in the evalu- 
ation of sub-leading pieces of two-nucleon effects in the 
process 7d^ 7r°d. Since u and w both display the same 
short-distance behavior, we arrive at the same conclusion 
as for the matrix elements considered above: for n = 1 
and n = 2 the result is convergent (in spite of oscilla- 
tions), see Figs. 01 and [SJ while for n > 3 it is divergent. 
For n = 1 and 2 evaluation with the wave functions of 
Ref. [13 gives: 



dr 



u(r)w(r) 



, u(r)w(r) „ 9 

dr v V = 0.156 fm" 2 



0.141 fm" 1 , (30) 
(31) 



Both of these numbers are consistent with the trend 
of the results for finite A, as shown in Figs. 01 and 
03 However, the second result will not be attained 
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FIG. 4: Cutoff-dependence of the integral J dr UA(r)wA(r)/r, 
where ma (wa) is the 3 Si ( 3 Di) radial deuteron wave function 
found for cutoff A in Ref. [l6|. The cross indicates the result 
of evaluation with the wave functions of Ref. 0, Eq. (EOJ. 



until A 20 fm , since approximately 20% of the 
final number accrues in the region between r — and 
r = 0.1 fm. 

Conclusion: In this paper we have shown explicitly 
that deuteron matrix elements of two-nucleon operators 
which are proportional to 1/r™ at short distances 
converge for n < 2 and diverge for n > 3 when they 
are evaluated using the leading-order deuteron wave 
function. We have given numerical evidence by explicit 
calculations using a sequence of leading-order deuteron 
wave functions corresponding to different ultraviolet 
cutoffs, and the results obtained in this way agree with 
an analysis based on the short-distance behavior of the 
deuteron wave function. 

As higher orders are calculated in the chiral series for 
the operator O the divergent case n = 3 will be reached 
(e.g. 7 d-+ 7r°d at 0(P 4 ) in X PT 0). Suppose this 
occurs at some chiral order m. If both the wave func- 
tion and the operator O are written as a chiral series, 
then this implies that an mth-order piece of the matrix 
element specifically the piece (ipi,o\0^ m '\ipi l o)i is di- 
vergent. This apparently mandates the presence of a con- 
tact interaction involving both nucleons and the external 
probe in 0^ m \ so that this divergence can be absorbed. 
However, in explicit calculations of a particular process 
to order m in the chiral expansion this divergence may 
cancel with contributions to J2J) due to corrections to 
the leading-order wave function. Such corrections will 
come from higher-order pieces in the chiral expansion of 
the NN potential, and will generate contributions to J2J, 
e. g., of the type S (V> (m) |O lo |VW • Therefore, the ap- 
pearance of a matrix element of an operator 1/r™ with 
ri > 3 at order m does not immediately indicate that a 
counterterm is needed at that order. However, if a coun- 
terterm is not permitted at that order by chiral symme- 
try and/or electromagnetic gauge invariance, then the 



0.3- 

A 

^ 0.1- 



10 

A[fm" 



15 



20 



FIG. 5: As for Fig. 0] but for / dr UA(r)wA{r)/r 2 . In this 
instance the cross is the result of Eq. 13H . 



various divergences at order m must cancel each other. 
This represents a constraint on the sum of all mechanisms 
that contribute to at that order. (For an example of 
such a cancellation in the context of pion production, see 
Refs. 0113.) 

Lastly, we note that if an NN potential that is 
more singular than one-pion exchange is iterated to all 
orders using the Schrodinger equation, as is done in 
Refs. 0, E1I3' then the behavior of the wave functions 
u and w will not follow the form 117|) . Indeed, the 
more singular the potential, the more convergent u and 
w will become 26]. Therefore, the conclusion of the 
previous paragraph is limited to an approach to higher- 
order calculations in the chiral effective theory where 
all corrections to the leading-order \ip) are evaluated 
in perturbation theory, and so the form (|17(l still rep- 
resents the dominant short-distance behavior of u and w. 
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